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ABSTRACT 

Bolocam began collecting science data in 2003 as the long-wavelength imaging camera at the Caltech 
Submillimeter Observatory. The planets, along with a handful of secondary calibrators, have been 
used to determine the flux calibration for all of the data collected with Bolocam. Uranus and Neptune 
stand out as the only two planets that are bright enough to be seen with high signal-to-noise in 
short integrations without saturating the standard Bolocam readout electronics. By analyzing all of 
the 143 GHz observations made with Bolocam between 2003 and 2010, we find that the brightness 
ratio of Uranus to Neptune is 1.027 ± 0.006, with no evidence for any variations over that period. 
Including previously published results at ~ 150 GHz, we find a brightness ratio of 1.029 ± 0.006 with 
no evidence for time variability over the period 1983-2010. Additionally, we find no evidence for 
time-variability in the brightness ratio of either Uranus or Neptune to the ultracompact HII region 
G34.3 or the protostellar source NGC 2071IR. Using recently published WMAP results we constrain 
the absolute 143 GHz brightness of both Uranus and Neptune to ~ 3%. Finally, we present ~ 3% 
absolute 143 GHz peak flux density values for the ultracompact HII regions G34.3 and K3-50A and 
the protostellar source NGC 2071IR. 

Subject headings: instrumentation: photometers, ISM: individual (G34.3, K3-50A, NGC 2071IR), 
planets and satellites: individual (Neptune, Uranus) 



1. INTRODUCTION 

The past two decades have seen vast improve- 
ments in the sensitivity and detector counts of broad- 
band (sub) millimeter imaging cameras, from single 
pixel rece ivers to backgroun d- limited kilopi xel ar- 
rays (e.g., iDuncan et all (I1990T); Holland et all (119991); 
iRunvan et all (120031); iDowell et all (120031); lHaig et all 
(12004: iCarlstrom et all (|2009[) : iSwetz et all (j2011l ); 
iHolland et all (|2006l )). Many of these instruments have 
relied on the planets, along with a handful of secondary 
calibrator s to o btain < 10% flux calibrations (e.g., 
IHill et all (120091) . IRudv et all (fl987[) IGriffin and Ortonl 
(|1993. hearafter G93f ). and iSandelf (|1994D 1 Alterna- 
tively, some of the large-scale survey instruments have 
used cosmic microwave background (CMB) measure- 
ment s to calibrate their da ta to even better precision 
(e.g..lRcichard t~ all (|2009f l) . However, calibrating from 
the CMB is only practical for instruments that image 
> 100 deg 2 , and therefore cannot be done in many cases. 
When planets and secondary calibrators are used, the 
accuracy of the flux calibration is often limited by uncer- 
tainties in the brightness of these c alibrators rather tha n 
by measurement uncertainties (e.g. JSavers et al.l ([2009)1 
Consequently, accurate measurements of the absolute 
brightness of the planets and secondary calibrators, along 
with a detailed understanding of any temporal variabil- 
ity in these brightnesses, is critical to maximizing the 
scientific potential of current and future (sub) millimeter 
cameras. 

Until recently, systematic uncertainties in the bright- 
ness of Uranus and Neptune have been ~ 5 — 10%, 
owing to uncertainties in the Martian brightness used 
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to deter mine the abso l ute brightness of e ach of these 
planets (iWrightJ (119761 ) . iRudv et all (|1987l l lOrton et all 
(|1986. hearafter 0861 ). G93). However, recently pub- 
lished WMAP results have significantly improved the 
constraints on the brightness of Mars (~ 1%) and 
Uranus (~ 3%) by calibratin g them relative to the CMB 
(jPage et al.l[2003t IliTetll] 127)091 : IWeiland et al.l[20Tll) . 
WMAP has also published CMB-calibra ted brightness 
meas urements of Neptune accurate to 8% (Weiland et al. 

EoUJ). 

Additionally, recently published results have shown 
evidence for a ~ 10% decrease in the bright- 
ness of Uranus at 8.6 and 90 GHz over the ~ 
20 ye ar period from the mid 1980s to the mid 
2000s (iKlein and Hofstadterl (12006. hearafter K106D and 
Kramer et all (|2008 hearafter Kr08l )l which have been 
attributed to changes in the relative viewing angle of the 
planet from earth due to its large obliquity of 82 deg. 
Note that, over the period of observations in K106 and 
Kr08 the SEP latitude of Uranus ranged from a min- 
imum of —82 deg in 1985 to a maximum of —7 deg 
in 2005. These data indicate that the brightness of 
Uranus increases as the magnitude of the sub-earth point 
(SEP) latitude increases (i.e., the south pole of Uranus is 
brighter than the equatorial regions). Further evidence 
for this scenario comes from resolved cm-wave images of 
Uranus, which indicat e that the poles ar e brighter than 
the equatorial regions (jOrton et alll2007D . 

2. OBSERVATIONS 

Bolocam is a large-format bolomctric camera with an 
8 arcmin field of view capable of observing at either 143 
or 268 GHz. Bolocam began collecting scientific data 
in early 2003 at the Caltech Submillimeter Observatory 
(CSO), and wa s commissioned as a facility instrument 
later that year (|Haig et al. 2004). In this manuscript we 
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TABLE 1 

Summary of 143 GHz Bolocam observations 



TABLE 2 

Bolocam 143 GHz brightness ratios 



date 


Uranus 


Neptune 


NGC 2071IR 


G34.3 


K3-50A 


2003/11 


6 


5 


9 








2004/11 


11 


20 











2006/11 


6 


3 


15 








2008/06 


2 


3 





5 





2008/07 


1 


2 





10 


2 


2008/11 


7 





8 








2009/10 


5 


21 


9 








2010/02 








18 


12 





2010/10 


6 


4 












date 


B Uranus 


B Uranus 


B N eptune 


B N eptune 


B «GC/G34 


B NGC/G3i 



Note. — The number of integrations for each flux calibra- 
tor for nine 143 GHz Bolocam observing runs between 2003 and 
2010. Note that there were five additional 143 GHz observing 
runs during that period, but those five observing runs included 
observations of only one flux calibration source, and are therefore 
excluded from our analysis. Note that the typical S/N of each 
integration is > 100. 

focus exclusively on 143 GHz Bolocam data, which was 
collected between 2003 and 2010 during nine separate 
observing runs (see Table [lj . All of the data collected 
with Bolocam has been calibrated using the planets (gen- 
er ally Uranus an d Neptune), along with sources given 
in Sandcll (1994). The Boloc am flux calibration p roce- 
dure is descr i bed i n detail in lLaurent et al.l (|2005[ ) and 
I Savers et al.l (|2009l ). and we briefly describe the details 
below. 

Fundamental to our flux calibration technique is the 
fact that Bolocam continuously monitors the operating 
resistance of the bolometers via the bias carrier ampli- 
tude. As the transmission of the atmosphere increases 
the optical load on the bolometers decreases and the op- 
erating resistance increases. Additionally, the responsiv- 
ity (in nV/Jy) of the bolometers increases as the bolome- 
ter resistance increases. Consequently, we are able to si- 
multaneously account for changes in atmospheric trans- 
mission and detector responsivity by fitting the flux cal- 
ibration as a function of bias carrier amplitude. In prac- 
tice, the technique works as follows. For a given observ- 
ing run, we generally observe multiple sources that have 
constant brightnesses. We then simultaneously fit the 
data for all of these sources according to 

Vi(Rj) = BiClijiot! + a 2 Rj), 

where Vi is the bolometer response to source i (in nV), 
Bi is the brightness of source i, fijj is the solid angle of 
source i during observation j, Rj is the bolometer resis- 
tance during observation j as measured by the bias car- 
rier amplitude, and a\ and «2 describe how the bolome- 
ter response changes as a function of Rj . 

Implicit in the above formula is the assumption that 
changes in the overall opacity of the atmosphere are com- 
pletely accounted for via a,\ and a2, which implies that 
the atmospheric transmission varies in a such way that 
the shape of the Bolocam bandpass through the atmo- 
sphere is constant. However, over the range of conditions 
where these data were collected (column depths of pre- 
cipitable water between 0.5 and 3.0 mm), variations in 
the atmospheric opacity do cause th e Boloc am bandpass 
to vary slightly (jPardo et al.ll2001aHb1 . [20051) ). For a typ- 
ical observing run, these bandpass variations in the at- 
mosphere add an additional uncertainty of ~ 0.2% to our 
measured brightness ratios, which is negligible compared 



2003/11 

2004/11 
2006/11 
2008/06 
2008/07 
2008/11 
2009/10 
2010/02 
2010/10 
total 



039 ± 0.045 
927 ± 0.065 
024 ± 0.010 
007 ± 0.022 
028 ±0.017 



1.029 ±0.026 1.016 ±0.042 



1.037 ±0.011 



1.006 ±0.009 
1.010 ±0.024 
1.038 ±0.030 
0.979 ±0.012 
0.995 ±0.012 



1.007 ±0.010 
1.028 ±0.015 
1.035 ± 0.032 

0.985 ± 0.008 



054 ± 0.047 
027 ± 0.006 



1.000 ±0.006 1.000 ±0.006 



Note. — Brightness ratios of Uranus to Neptune, 
Uranus to the secondary calibrators NGC 2071IR and 
G34.3, and Neptune to the same secondary calibrators. 
Note that the brightness ratios of the planets to the sec- 
ondary calibrators are normalized to the mean brightness 
ratio since the solid angle of the secondary calibrators is 
not precisely known. The brightness ratio of NGC 2071IR 
to G34.3 is constrained using the 2010/02 data. 

to our measurement uncertainties of 1 — 5%. 

Note that we have assumed that the angular size of the 
source is much smaller than the angular size of the Bolo- 
cam point-spread function, which is true for all of the 
sourcefl Since the solid angle of both Uranus and Nep- 
tune varies with observation epoch, we compute the value 
of Ojj separately for each integration of these planets us- 
ing the James Clerk Maxwell Telescope (JCMT) Fluxes 
prograuQ. For the secondary calibrators, the value of 
is assumed to be constant. See Table [2] and Fig- 
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3. RESULTS 



Seven Bolocam observing runs contain observations of 
both Uranus and Neptune. In each run, we measure the 
brightness ratio of Uranus to Neptune with < 5% sta- 
tistical uncertainty. If we assume that the brightness 
ratio is constant over the period from 2003 to 2010, then 
we find that B Uranus / B Neptune = 1.027 ± 0.006 with 
X 2 /DOF = 4.5/6. A linear fit versus year to the data 
reduces the value of x 2 by 0.9, indicating that a linear fit 
is not required (a F-test shows that 39% of realizations 
of constant data would have yielded a larger F-ratio, im- 
plying that a constant fit to our data is adequate) . Com- 
paring to published results at similar frequencies, 086 
find a brightness ratio of 1.016 ± 0.038 using 150 GHz 
data collected between February 1983 and March 1984 
at the National Radio Astronomy Observatory 12-m an- 
tenna at Kitt Peak, and G93 find a brightness ratio of 
1.043 ±0.016 using 156 GHz data collected at the JCMT 
in May 1990/1992. The combined ~ 150 GHz data set of 
086, G93, and Bolocam is well described by a constant 
brightness ratio of B Uranus / B Neptune = 1.029 ± 0.006, 
with a x 2 /DOF = 5.4/8. In this case a linear fit versus 

3 The semi-diameters of the sources in arcsec are ~ 2 (Uranus), 
~ 1 (Neptune), ~ 6 (NGC 2071IR), ~ 3.5 (G34.3), and ~ 4 (K3- 
50A) , all of which are small compared to Bolo cam's 59 arcse c full- 
width at half-maximum point-spread function (Sandcll 1994]). The 
secondary calibrators are large enough that our point-like assump- 
tion will cause us to underestimate the total surface brightness of 
these objects by < 1% (see Section [3). Uranus and Neptune are 
small enough that the correction is negligible (< 0.1%). 

4 http:/ /www .jach.hawaii.edu/jac-bin/planetflux. pi 
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Fig. 1. — From left to right, Bolocam measurements of the brightness ratio of Uranus to Neptune, the relative brightness ratio of Uranus 
to NGC 2071IR or G34.3, and the relative brightness ratio of Neptune to NGC 2071IR or G34.3 for the period from 2003 to 2010. The 
red line shows the best-fit mean value of the brightness ratio (equal to 1 by definition for the center and right plots), and the grey band 
represents the la confidence region. Our data show no evidence for temporal variations in any of these brightness ratios. 



year reduces the value of x 2 by 0.1, and a F-test shows 
that 76% of realizations of constant data would have pro- 
duced a larger F-ratio. In addition, the linear fit provides 
a formal 95% confidence level upper limit of 5.3% on the 
magnitude of the variation over the 28-year period 1983 
to 2010 (SEP latitudes of -82 to +10). 

In addition to searching for variations in the brightness 
ratio of Uranus and Neptune, we also searched for tem- 
poral variations in the brightness of each planet relative 
to sources that are known to have const ant brightnesses, 
NGC 2071IR and G34.3 (jSandelll [l99l . Six observing 
runs contained observations of Uranus and either NGC 
2071IR or G34.3, and five observing runs contained ob- 
servations of Neptune and either NGC 2071IR or G34.3. 
Note that we constrained the brightness ratio of NGC 
2071IR and G34.3 based on observations of both ob- 
jects in February 2010. A constant fit of the bright- 
ness ratios yields a x 2 /DOF= 6.3/5 for Uranus and a 
X 2 /DOF= 9.1/4 for Neptune. We again calculate the F- 
ratio for a constant fit to the brightness ratios compared 
to a linear fit versus year and again find that there is 
not compelling evidence for a linear fit (22% of realiza- 
tions of constant data would have a larger F-ratio than 
the Uranus data and 21% of realizations of constant data 
would have a larger F-ratio than Neptune data). 

Since we are unable to measure any temporal varia- 
tions in the brightnesses of Uranus or Neptune, we use 
the WMAP measure ments of the brightn ess tempera- 
ture of each planet (jWciland c t al.1 120111 ). along with 
our measurement of the brightness ratio (1.027) to de- 
termine the 143 GHz brightness temperature of each 
planet. In G93 the authors present models for the tem- 
perature spectra of each planet based on observations 
between 90 and 850 GHz, and we have used these mod- 
els to extrapolate the 94 GHz WMAP measurements to 
143 GHz. G93 state that these models are accurate to 
< 2 K when absolute calibration uncertainties are not 
included, and we have therefore included a 2 K uncer- 
tainty in the model-based extrapolation of the absolute 
94 GHz WMAP measurement to 143 GHz. We find that 
the brightness temperature of Uranus is 106.6±3.5 K and 
the brightness temperature of Neptune is 103.8 ± 3.4 K, 
where the uncertainties include the WMAP and Bolo- 
cam measurement uncertainties, along with the quoted 
uncertainties in the G93 models. These brightness tem- 
peratures indicate that the G93 models need to be scaled 



TABLE 3 

Absolute 143 GHz Brightness Values 



object 


brightness 


Uranus 


106.6 ±3.5 K 


Neptune 


103.8 ±3.4 K 


NGC 2071IR 


1.315 ±0.043 Jy 


G34.3 


13.87 ±0.47 Jy 


K3-50A 


7.387 ±0.244 Jy 



Note. — Absolute 143 GHz brightness values 
for calibrators observed by Bolocam. Uranus and 
Neptune are given in surface brightness units (K) 
and the secondary calibrators are given in units of 
peak flux density per beam (Jy). 



by 0.931 ± 0.031 (Uranus) and 0.946 ± 0.031 (Neptune). 
Since the G93 brightness models of Uranus and Nep- 
tune were n ormalized to t he Ma rtian calibration model 
described in [Griffin et all (11986D (based on results from 
lWrightJ (ll976ira nd Ulich T(| 19811 )). our results indicate that 
the IGriffin et al.l (jl986[ ) model over-predicts the bright- 
ness of Mars by ~ 5 — 7%. These scaling s are consistent 



wit h the results o f lHalverson et all ((2009), who find that 
the lWrighl] (|1976| ) model should be scaled by a factor of 
0.922 ± 0.037 based o n the WMAP Ma rtian brightness 
measurements given in IHill et al.l ff2009h . 

Using these absolute brightness temperatures of 
Uranus and Neptune, we find peak flux densities per 
beam of 1.315 ±0.006 ±0.043 Jy (NGC 2071IR), 13.87± 
0.10 ± 0.46 Jy (G34.3), and 7.387 ± 0.009 ± 0.244 Jy 
(K3-50A) for the three secondary calibrators (see Ta- 
ble [3]) . The first value represents our measurement un- 
certainty on each peak flux density, while the second 
value represents our 3.3% uncertainty in the brightness 
temperatures of the planets. Given the source sizes in 
ISandelll (jl994[ ). along with Bolocam's 59 arcsec FWHM 
PSF, these values should be withi n 1% of t he to tal flux 
density for each source. Note that ISandelll (|1994l ) found 
the following peak flux densities per 27 arcsec beam for 
these sources at 150 GHz: 1.5 ± 0.2 Jy (NGC 2071IR), 
12.5 ± 0.4 Jy (G34.3), and 6.5 ± 0.2 Jy (K3-50A), where 
the uncertainties do not include systematic uncertainties 
in the planet models used to calibrate the data. Our 
measured peak flux densities are ~ 10% higher for both 
G34.3 and K3- 50A (plus an ad ditional factor of ~ 7% due 
to the fact that lSandell (|1994l) used the G93 model in or- 
der to determine his absolute calibration, which we have 
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Fig. 2. — Uranus brightness measurements relative to the average 
brightness of Uranus at a range of frequencies between 1980 and 
2010. The black points show ~ 150 GHz data from 086, G93, and 
Bolocam, the red triangles show 90 GHz data from Kr08, and the 
green squares show 8.6 GHz data from K106 (taken from Table 2 
of each manuscript and rebinned for visualization). 

found to overestimate the brightness of Uranus by ~ 7%) . 
It is not clear why our measured peak flux densities are 
higher for G34.3 and K3-50A, although the known ex- 
tend ed emiss i on in both sources, coupled with the fact 
that iSandelll ([1994) made single-pixel chopped photom- 
etry measurements at 27 arcsec resolution, may be the 
cause. 

4. DISCUSSION 

We find that Uranus and Neptune behave as ideal 
sources for flux calibration at 143 GHz, with no evidence 
for temporal brightness variations. For Uranus, these re- 
sults are in contrast to the lower frequency measurements 
of K106 and Kr08, who find ~ 0.5 percent /year varia- 
tions in the brightness temperature of Uranus at 8.6 and 
90 GHz. Our data, combined with the measurements of 
086 and G93, place a 95% confidence level upper limit of 
0.19 percent/year on the magnitude of variations in the 
brightness temperature of Uranus at ~ 150 GHz over the 
same period. See Figure [5J A physical interpretation of 
the temporal variations in the brightness of Uranus seen 
at lower frequencies by K106 and Kr08, in combination 
with our static 143 GHz results, is beyond the scope of 
this manuscript, which is intended to quantify the mag- 
nitude and stability of the brightness of Uranus for the 



purposes of using it as a calibrator at 143 GHz. However, 
we do note that the combined results are not necessarily 
inconsistent, given that higher frequency observations of 
Uranus probe higher altitudes in the atmosphere (Kr08). 

5. CONCLUSIONS 

Using Bolocam data collected between 2003 and 2010 
we have tightly constrained the 143 GHz brightness ratio 
of Uranus and Neptune (1.027 ± 0.006), and we find no 
evidence for temporal variations in the 143 GHz bright- 
ness temperature of either planet over that period. Com- 
bining our results with those of 086 and G93, we find 
no evidence for 143 GHz brightness variations in either 
planet over the period from 1983 — 2010, and place a 95% 
confidence level upper limit of 5.3% on the magnitude of 
brightness variations over the 28 year period from 1983 to 
20 10. By extrapolating the WMAP 94 GHz results given 
in iWeiland et al.l (|2011l) to our observing band using the 
brightness models presented in G93, we are able to con- 
strain the absolute 143 GHz brightness temperature of 
each planet to ~ 3%. Additionally, we determine ~ 3% 
absolute 143 GHz peak flux densities for the ultracom- 
pact HII regions G34.3 and K3-50A and the protostellar 
source NGC 2071IR. 
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